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Abstract 

In the framework of the relativistic quasipotential quark model the mass spec- 
trum of baryons with two heavy quarks is calculated. The quasipotentials for in- 
teractions of two quarks and of a quark with a scalar and axial vector diquark are 
evaluated. The bound state masses of baryons with = ^ , | are computed. 

1 Introduction 

The investigation of properties of hadrons, containing heavy quarks, is of a great interest 
for understanding the dynamics of quark and gluon interaction. To calculate the meson 
and baryon mass spectra QCD sum rules potential models ||, |], lattice QCD 

heavy quark effective theory (HQET) 0, §], the method of vacuum correlators p 
and some other approaches are widely used. Presently at the LHC, 5-factories and 
the Tevatron with high luminosity, several experiments have been proposed, in which a 
detailed study of baryons containing two heavy quarks can be performed. The possible 
quark composition of such baryons looks as follows: (ccg, cbq, bbq), where c and b are 
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heavy quarks (Q) and q denotes a light u, d and s quark. In this connection theoretical 
predictions for heavy baryon masses acquire important significance. The properties of 
heavy hadrons (containing h and/or c quarks) are essentially different from those of light 
hadrons composed from m, d and s quarks. This difference originates from the fact that the 
heavy quark mass strongly exceeds the scale of QCD interaction, Is.qcd ~ 300 — 400 MeV, 
namely mc,f, ^ ^qcd- Therefore, while investigating processes with heavy quarks one 
should expect that the limit mq —>■ oo within QCD yields good quantitative results, which 
may be further improved in HQET by taking into account l/mg corrections. Baryons, 
containing two heavy quarks (QQq), may be considered as a localized source (QQ) of 
the colour field, in which the light quark moves. Interaction forces between heavy quarks 
then lead to the formation of a two-particle bound state of a QQ-diquark, the scale of 
which is determined by the quantity l/mq which is small compared to the QCD scale 
1/Aqcd- Thereby it seems justified to treat the heavy diquark as a pointlike object with 
definite colour, spin and mass and to ignore the configuration with a gQ-diquark which 
size is of order l/rUg ~ 1/Aqcd- The baryon formation then occurs as a result of a 
diquark interaction with the light quark. Here in the framework of the quark-diquark 
approximation we calculate the mass spectrum of heavy baryons on the basis of a local 
Schrodinger-like quasipotential equation |[T0[] : 

where the relativistic reduced mass is 

E,E2 _M^-{ml-mlf 

M'^ — ml + m\ M^ — ml + ml 
E-i = , Eo = , E^ + Eo = M, 

and the center of mass system relative momentum squared on the mass shell reads 

.2,..^ _ [M' - (mi + m2)^][M^ - {m, - m^)^] 



with mi 2 the masses of the constituent particles. In refs. [|rT|, |1^] the operator of quark- 
antiquark interaction has been constructed and the mass spectra and decay rates of mesons 
have been investigated. The relativistic quasipotential quark model gives in the meson 
sector results which nicely agree with experimental data. Relativistic effects play an 
important role in describing properties of quark bound states and may be consistently 



taken into account within the quasipotential method 13 . 



2 Diquarks in the relativistic quasipotential model 

The concept of diquarks was introduced at the very beginning of the development of the 
quark model of hadrons in order to describe baryon properties [^. Diquarks represent 
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by themselves two-particle clusters, which are produced inside the three-quark systems 
as a result of quark interactions in which the spin-dependent forces play an important 
role [|l4| , p!5[] . The attractive forces between two quarks act if the quarks are in the 



antisymmetric colour state. The interaction potential at large distances in the quark- 
quark system is similar to that in the quark-antiquark system since the baryon size is 
practically the same as the meson size. In the ground state the diquarks are two-particle 
bound states of quarks in an antisymmetric colour state with zero angular momentum 
and definite flavour and spin. In the case of indentical quarks the diquark has spin S = 1, 
whereas for quarks of different flavours the state with spin S' = is also possible. 

Now we construct the interaction quasipotential for two heavy quarks Q, which enters 
in eq. (1). The interaction operator in the QQ system contains a perturbative part, which 
is defined at small distances in the quasipotential approach by the one-gluon exchange 
amplitude and can be evaluated within QCD. Moreover there is a nonperturbative term 
in the operator V{p, q, M), which cannot be obtained consistently within QCD and leads 
to the quark confinement at large distances. 

The quasipotential V{p, q, M) of eq. (1) is constructed with the help of the off-mass- 
shell scattering amplitude, projected onto the positive energy states. Then the perturba- 
tive part of the quasipotential can be represented in the form: 



VqQiP^ ^ = Mp)M-P) I ^ D^^{k)-i'(Y2Ui{q)u2{-q), (4) 

4Vei(p)ei(g)e2(p)e2(gj 

where D^j_y{k) is the gluon propagator conveniently taken in the Coulomb gauge 

D°°(fc) = -1, D^^(A:) = -1(^5--^), D'^ = D^ = Q, (5) 

k = p — q, and u{p) is the Dirac spinor, 

= (6) 

\s{'p)+m/ 

In this case the colour factor, which includes the colour part of the baryon wave function 
£ijk/ k = 1,2, 3) has the value: 



^j^^ijk^nk ^j^^lmnTijT^j — ^, (7) 



where T" are the generators of the colour S'f/c(3) group. It means, that in the two- 
quark system the exchange of a nonabelian gauge boson results in an attractive Coulomb 
potential at small distances which differs only by a factor 1/2 from the similar potential 
in the quark-antiquark system. Substituting the spinor (6) into eq.(4) and carrying out 
some manipulations with spinor convolutions, we obtain, in the same way as in WM, the 



one-gluon exchange quasipotential including relativistic corrections. Since we consider the 
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ground state of diquarks, the angular momentum of the QQ system is zero. As it has been 
remarked in ref . []T6| , the quark spin interaction can result in the existence of short-range 



correlations in the quark-quark subsystem. Therefore, neglecting corrections 0{p'^ /rn^) 
we can represent the perturbative part of the quasipotential in the configuration space in 
the form: 

= -\- + ^^{a,a,)5{f), (8) 
6 r 9mim2 

where ^^(r) is the running coupling constant and ai^2 are the quark spin matrices. The 
colour hyperfine interaction is attractive if the total spin of two quarks S = 0, and 
repulsive in the case S = 1. 

In the approximation of one-gluon exchange the quark-quark potential in a baryon 
is equal to one half of the quark-antiquark potential in a meson. Calculations based on 
the baryon Wilson loop area law in QCD show that this relation holds even true for 



nonperturbative interactions. Thus, while constructing the long-range nonperturbative 
part of the quark interaction operator we used the known "1/2 rule", i. e. Vgg = \Vqq p 
and chose this part of the quasipotential in the standard linearly rising in r form. As a 
result in this approximation the complete interaction operator in the quark-quark system 
is chosen in the form: 

VQQ{r) = + \{Ar + B) + ^^{a,a,)6{r), (9) 

6 r I 9mim2 



where the parameters A = 0.18 GeV^ and B = —0.3 GeV had been fixed previously [11 



T2| in calculating the meson mass spectrum. The numerical solution of the quasipotential 
equation (1) with interaction operator (9) yields the mass spectrum of the s-wave axial 
vector and scalar diquarks. The results are presented Qin Table 1. Here the usual values 
of quark masses (the same as for mesons) were used ||Tl|, |12[ : 

m„ = mrf = 0.33 GeV, m, = 0.5 GeV, = 1.55 GeV, nib = 4.88 GeV. 



Table 1: Masses of s-wave scalar and axial vector diquarks 



Diquark 


qq 


qs 


qc 


qb 


ss 


sc 


sb 


cc 


cb 


bb 


M(GeV) 
(scalar diquark) 


0.85 


1.04 


2.08 


5.39 




2.27 


5.57 




6.52 




M(GeV) 
(axial vector diquark) 


1.02 


1.15 


2.12 


5.41 


1.32 


2.29 


5.58 


3.26 


6.52 


9.79 



The masses of light diquarks are presented for completeness where q = u,d. 



4 



3 Quasipotential for the interaction of a quark with 
a scalar diquark 

Now to calculate the baryon mass spectrum in the quark-diquark approximation it is 
necessary to construct the quark-diquark interaction operator. In the case of a scalar 
diquark the perturbative part of the quasipotential is described by the Feynman diagram 
in Fig. 1. 



Pi = {Ei,p) 



P2 = {E2, -p) 




qi = {Ei,q) 



q2 = {E2, -q) 



Figure 1: Feynman diagram of the one-gluon interaction in the quark-diquark system. 
The corresponding expression for the interaction operator has the form: 

V^TsniP, a = , ,A P2 + Q2rD,.{k)u{q)Yu{p), (10) 

4vei(pjei(g)e2(pje2(gj 

where the interaction vertex of a scalar diquark with a gluon is given by the factor 
igsT^j{p2 + 12)^- Using the standard representation of the Dirac spinors (6), we get for 
^q+SD up to terms 0{p^ /w?) the following expression: 



V^;%{p.q)--^' 



1 + 



— b"^ if — b"^ ff — b"^ , 

+ , 2 + , 2 +- (11 

4mf mim2 



E1E2 

-(P^) 17-^ + 7-^ + ——) + iSi\pxq\-^(— + —]\. 

\Ami mim2 J 4mim2 2mi\mi fn,2j) 

Then we add to this expression in the configuration space the confining potential 

C5r*-(r) = Ar + S, (12) 

which is the same as in the quark-antiquark system. Thus the complete quasipotential, 
which is used to calculate the mass spectrum of baryons with = |^ can be represented 
as follows: 
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2iial (I 1 4 \ 4a,, / 1 1 1 



3r^ \^m^ ' m2 mim2 J 3r^ ^ ' ' \^4mf 4m2 mim2 

where m2 is the scalar diquark mass, mi is the hght quark mass and fi is their reduced 
mass. The numerical results for the mass spectrum of the ground states of baryons with 
= 1^ on the basis of eqs. (1) and (13) are given in Table 2. 

4 Quasipotential for the interaction of a quark with 
an axial vector diquark 

As in the previous section, the perturbative part of the interaction operator for spin 1/2 
and spin 1 particles is given by the diagram in Fig. 1. It includes the matrix element of 
the chromomagnetic current operator between states of the axial vector diquark. This 
matrix element is defined by three chromomagnetic form factors. Neglecting the chromo- 
quadrupole moment of the axial vector particle, we can represent the four vector of the 
chromomagnetic current in the form |]TB|: 

fip, q) = tg^C{Gi{k'){p2 + q2Y + ^^[W>'{mS2^2) - {S2^2)W^{n)]]i, (14) 

A2 = -q-\ -C/2 



1712 \ E2+m2j 

where 5*2 is the spin operator of a vector particle, IV^ is the four vector of the relativistic 
spin (the Pauli-Lubansky vector) with the components: 

W%p2) = -S2P, W{p2)=m2S2 + ^^^, (15) 

E2 + m2 

and ^ is the axial vector particle wave function (polarization vector) in its rest frame. 
Then the one-gluon exchange quasipotential takes the form: 

v',:iu{p. ^ = -^==l^—=D,.{k)miMp} (16) 

4\/ei(p)ei(gje2(p)e2(gj 



Giie){p2 + q2Y + G2{e)[—W^'{p2){S2^2) - — (5^2A2)W^'^(P2)] 

m2 m2 

Now we expand eq. (16) with an accuracy up to terms 0{jP' /m^). The zeroth and spacial 
vector components of the current (14) contain the following terms: 

— Vr°(p2)(5^2A2) - — (52A2)iy°(p2) = — 5^2[P X g], (17) 

1712 ^2 ^2 

— W{p2){S2^2) - —{S2^2)W{p2) = l[{p- q) X ^2], (18) 

m2 Tn2 
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which produce fine and hyperfine sphttings of energy levels. The chromomagnetic form 
factor of the axial vector particle G2{k'^) at /c^ = defines its total chromomagnetic 
moment: (^2(0) = /iy = 1 + C- Setting here the anomalous chromomagnetic moment of 
the axial vector diquark C = 1 ^-^id 6*1(0) = 1, we obtain the following expression for 



Tr-pert. 
^q+AD- 



q+AD[P,<l)- 3 ^2 \['^E,E2 Ami ^rnl ^ m^m^ 



.ir^ + ir^ + + -^ + '^^Ap X — — + — ] 

4m| mim2 Amim2 2,mi \mi m2 J 



+tS2\px qIt^I— + —]- -^[iS,S2)k' - iS,k)iS2k)]]. 

21712 \m2 ^1 / mim2 J 

The nonperturbative part of the interaction potential for a quark and an axial vector 
diquark is taken in the form (12). Transforming (19) into configuration space, we get the 
complete expression of the interaction operator for the system q + AD: 

Aas A b"^ \ Anas 4a<j f 1712 \ LSi 



2ua2 f 1 1 4 \ 4a,, / I 1 1 \ 



3^2 ^2 777,-^777,2 y 3r^ y4mf 41712 rnim2^ 

6mim2 \2m2 / r'^ 9mim2 6mim2r-^ 

The potential (20) has been used for calculating the mass spectrum of baryons with = 
I ,| on the basis of the local quasipotential equation (1). Corresponding numerical 
results are presented in Table 2. 



5 Conclusions 



In this paper we have carried out the calculation of the mass spectrum of baryons with two 
heavy quarks in the quark-diquark approximation on the basis of a local quasipotential 
equation. As it follows from Table 2, our results are in good agreement with calculations 
made within other approaches. One can see from Table 2 that our predictions are es- 
pecially close to those obtained in |T^. The approach used in [|r^, ^ is based on some 
semiempirical regularities of the baryon mass spectrum following from experimental data 



and the constituent quark model. The authors of [|T^, |20[ found a successful parametriza- 
tion of the interpolating curves which describe the known meson and baryon masses, and 
give predictions for some yet unobserved mesons and baryons. Our results for the Be me- 



son mass spectrum |12| are also rather close to those of |20|. Our only remark is that one 
should take into account the mixing of ^Pi and ^Pi Be states. The quoted authors use in 
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Table 2: Masses of baryons (in GeV) containing two heavy quarks. [QQ] denotes the 
diquark subsystem with the antisymmetric spin wave function, {QQ} denotes the diquark 
subsystem with the symmetric spin wave function (g = u,d). 



Notation 


Quark 




Mb 


Mb 


Mb 




content 




(our results) 


[in 

[191 




'CC 


{ccjq 


1 + 


3. 00 


3.66 


3.61 


^cc 


{cc\q 


2 


O O 1 

o.ol 


3.74 


3.68 


^ ^cc 


{cc}s 


1 + 
2 


3.7o 


3.74 


3.71 




{cc}s 


3 + 


o on 
o.oy 


3.o2 


3. lb 


^bb 


{bb}q 




10.23 


10.34 




^bb 


{bb}q 


2 


10.28 


10.37 




^bb 


{bb}s 


1 + 
2 


10.32 


10.37 




0* 
"fob 


{bb}s 


3 + 
2 


10.36 


10.40 




^cb 


{cb}q 


1 + 
2 


6.95 


7.04 




^'cb 


[cb]q 


1 + 


7.00 


6.99 






{cb}q 


2 


7.02 


7.06 




^cb 


{cb}s 


1 + 
2 


7.05 


7.09 




^'cb 


[cb]s 


1 + 


7.09 


7.06 




^Ib 


{cb}s 


2 


7.11 


7.12 





particular the Feynman-Hellmann theorem in order to obtain the monotonic decrease of 
the hadron mass spectrum with respect to quark masses. Here they use an oversimplified 
assumption about the interquark potential, namely the complete neglection of the spin- 
independent part of the Breit-Fermi potential as being unimportant. This assumption 
leads finally to some inequalities between quark masses and meson masses. These in- 
equalities are violated in many constituent quark models (including ours) in which widely 
adopted quark masses are used. From our point of view there is nothing serious in these 
contradictions and the recovery of the appropriate spin-independent (flavour-dependent) 
part of the quark potential |]T2[ will eliminate them. Thus we may conclude that our 
model correctly describes some important features of the baryon mass spectrum. 

There exists yet another type of baryons with one heavy and two light quarks for 
which the quark-diquark approximation gives reasonable results. Thus, in Ref. Qqq 



baryons were treated by considering the interaction of the diquark system with the third 



quark on the basis of additional quark exchange forces p2|, Bl, E4| . The resulting effective 



lagrangian, which incorporates heavy quark and chiral symmetry, describes interactions 
of heavy baryons with Goldstone bosons in the low energy region. As an interesting 
application, the Isgur-Wise form factors have been calculated. 
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